Orbital selective and tunable Kondo effect of magnetic adatoms on graphene: 
Correlated electronic structure calculations 
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We have studied the effect of dynamical correlations on the electronic structure of single Co 
adatoms on graphene monolayers with a recently developed novel method for nanoscopic materials 
that combines density functional calculations with a fully dynamical treatment of the strongly 
interacting 3d-electrons. The coupling of the Co 3d-shell to the graphene substrate and hence the 
dynamic correlations are strongly dependent on the orbital symmetry and the system parameters 
(temperature, distance of the adatom from the graphene sheet, gate voltage). When the Kondo effect 
takes place, we find that the dynamical correlations give rise to strongly temperature-dependent 
peaks in the Co 3d-spectra near the Fermi level. Moreover, we find that the Kondo effect can be 
tuned by the application of a gate voltage. It turns out that the position of the Kondo peaks is 
pinned to the Dirac points of graphene rather than to the chemical potential. 

PACS numbers: 



I. INTRODUCTION 

Since its recent discovery, Graphene — a monolayer of 
graphite — has become the subject of intense research 
due to its peculiar electronic properties^— . The unusual 
properties of graphene are largely a result of the linear 
dispersion of electron bands at low energies and its two- 
dimensionality. This allows to describe the charge car- 
riers in graphene as a two-dimensional gas of massless 
Dirac-fermions. The unusual electronic properties and 
the fact that its chemical potential can easily be tuned 
by a gate voltage make graphene an ideal experimental 
testing ground for exotic physics as well as a promis- 
ing basis for novel nano-electronics devices. For recent 
reviews of this rapidly developing field and additional 
references, see e.g. Refs. 

Graphene is also a good testing ground for exotic 
Kondo physics^: First, the linearly vanishing DOS of 
graphene favours the formation of a local moment 9 which 
is a necessary condition for Kondo physics. Also vacan- 
cies and edge defects in graphene have been predicted to 
give rise to magnetic moment a 10 ' 11 . Very recently, the 
Kondo effect has actually been observed at local defects 
in graphene providing evidence that defects in graphene 
are indeed magnetic^. On the other hand, as a result 
of the linearly vanishing DOS, a magnetic moment cou- 
pled to a graphene layer can only display a Kondo ef- 
fect if the strength of the coupling is strong enough^— . 
Furthermore because of the existence of two inequivalent 
Dirac cones, it has multiple channels which can in prin- 
ciple screen the impurity spin and could therefore lead to 
a multichannel Kondo effect and hence non-fermi liquid 
behaviour— ~—. Kondo effect and the possibilty of mul- 
tichannel Kondo effect in graphene have been discussed 
in a number of recent publications 2 ^ - — . 



Moreover, it is possible to probe the Kondo effect of 
magnetic adatoms on graphene by scanning tunneling mi- 
croscopy (STM). Similar to STM spectroscopy of mag- 
netic adatoms and molecules on metallic surfaces 2 ^—, 
the Kondo effect could show up as Fano lineshapes in 
the conductance vs. bias voltage characteristics of the 
STM or as additional temperature dependent peaks in 
the STM spectra. But in comparison with normal metal- 
lic surfaces the Fano lineshapes for graphene are expected 
to be more asymmetric and strongly dependend on the 
symmetry of the d-orbitals of magnetic atom 2 ^. 

The goal of this work is to perform a realistic study 
of the Kondo effect for a single Co adatom on an other- 
wise perfect graphene sheet by using a recently developed 
approach for nanoscopic systems that combines ab initio 
electronic structure calculations on the level of density 
functional theory (DFT) calculations with a full dynam- 
ical treatment of the strongly interacting 3d-electrons2= 
A realistic treatment allows us to investigate under which 
conditions the Kondo effect takes place in a given orbital 
and the detailed spectral functions that arise from the 
Kondo effect in a semi-metallic substrate. 

This paper is organized as follows: In Sec. |TI]we briefly 
describe our recently developed GGA+OCA method 31 
for the quasi ab initio description of the electronic struc- 
ture of nanoscopic systems fully taking into account 
dynamical correlations by combining DFT calculations 
with a dynamical treatment of the strongly interacting 
3(i-electrons in the so-called one-crossing approximation 
(OCA). 

In Sec. IIIII we present DFT electronic structure calcu- 
lations in the generalized gradient approximation (GGA) 
for determining the most favourable adsorption site for 
the Co atom on the graphene sheet. In agreement with 
previous results by other group o 24 ' 25 ' 39 we find that the 
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hollow site is the most favourable adsorption site. Us- 
ing hybrid functional calculations which generally predict 
better geometries for molecules we have then optimized 
the adsorption height of the Co atom at the hollow site 
of the graphene. Using GGA we have also calculated 
the charge and spin of the Co atom in dependence of 
the distance to the graphene sheet for the three differ- 
ent adsorption sites. The plots indicate that depending 
on the distance the Co atom can be either in the local 
moment regime or in weak correlation regime. Since the 
potential energy curves are very shallow it seems feasible 
that the correlation regime be controlled via the distance 
of the Co atom to the graphene substrate in an actual 
experiment. 

In Sec. IIVI we study the impact of dynamical cor- 
relations of the Co 3d-electrons on the elctronic struc- 
ture of the Co atom and the graphene sheet using the 
GGA+OCA method for nanoscopic systems as described 
in Sec. [TTJ We find that the dynamical correlations of 
the Co 3d-electrons give rise to strongly temperature- 
dependent peaks near the Fermi level in the spectral den- 
sity of the Co 3d-shell but also in the spectral density of 
the p z -orbitals of nearby carbon atoms. Furthermore, we 
find that the Co 3d-spectra depend quite strongly on the 
exact energies of the Co 3<i-levels. Accordingly, the ap- 
plication of a gate voltage which changes the chemical 
potential of the graphene substrate has a similar effect 
on the spectra. Hence the dynamic correlations and the 
Kondo effect can be controlled by the application of a 
gate voltage. 

Finally, in Sec. |V]we conclude this work with a general 
discussion of the results in the light of Kondo physics. 



II. GGA+OCA METHOD 

We consider a single Co adatom adsorbed at the hollow 
site of an otherwise perfect graphene sheet as shown in 
Fig. [3] To proceed we divide the system into two parts. 
The central region C which contains the Co adatom and 
a number of carbon atoms is embeded into the host ma- 
terial H given by a perfect graphene sheet. The central 
region C has to be chosen big enough so that the elec- 
tronic structure outside of the region can be assumed to 
be that of a perfect graphene sheet. 

As a first step a mean-field description of the electronic 
structure of the system is obtained from DFT calcula- 
tions on the level of GGA using the CRYSTAL06 ab initio 
electronic structure codeSZ- together with the 6-31G gaus- 
sian basis set: By setting up a 2D lattice of supercells of 
the C region we obtain an effective one-body Hamiltonian 
He for the C region from the Kohn-Sham Hamiltonian of 
the converged GGA calculation. Analogously, we obtain 
an effective mean-field description of the host material 
from a GGA calculation for a supercell of clean graphene 
corresponding to the C region. 

The Green's function of the central unit cell C is then 



given by: 

GcM = (w + /i - H c + h dcc - E M (w) - ShM)- 1 (1) 

Sh(uj) is the embedding self-energy which describes the 
dynamic hybridization between the central region C and 
the host material H. The embedding self-energy is ob- 
tained from the GGA calculation of the perfect graphene 
sheet as explained in App. [X] 

E3d(w) is the electronic self-energy describing the 
dynamic correlations of the strongly interacting 3d- 
electrons. The strong interactions of the Co 3d-shell 
are captured by adding a Hubbard-like interaction term 
to the one-body Hamiltonian within the correlated 3d- 
subspace: 

HU = ^ ]C U » kl ^iajjvjl^dka! (2) 

cr l cr 2 

Uijki are the matrix elements of the effective Coulomb 
interaction of the 3d-electrons which is different from the 
bare Coulomb interaction due to the screening by the 
conduction electrons. Here we take a model interaction 
taking into account only the direct Coulomb repulsion U 
between electrons and a Hund's rule coupling J. For 3d 
transition metal elements in bulk materials the repulsion 
U is around 2-3 eV and J is around 0.9 eV— . Due to 
the lower coordination of the Co atom on the graphene 
sheet the screening should be reduced compared to its 
bulk value. Here we take {7 = 5 eV and J = 0.9 eV as in 
our previous work^i. 

The Coulomb interaction within the correlated 3d sub- 
space has already been taken into account on a static 
mean-field level in the effective KS Hamiltonian of the 
central unit cell C. Therefore the KS Hamiltonian within 
the correlated subspace hj£f has to be corrected by a 
double-counting correction term, i.e. = h^ d s — h^cc- 
Here we use the standard expression, 

h dcc = [U(N 3d - 1/2) - J(N 3d - l)/2] x l 3d (3) 

where is the identity matrix in the Co 3d-subspace, 
and N^d is the occupation of the impurity 3d-shell2£. 

In order to compute the electronic self-energy £ 3( j we 
have to solve the generalized Anderson impurity problem 
given by the strongly interacting Co 3d-electrons. The 
impurity problem is completely determined by the inter- 
acting Hamiltonian of the Co 3d-shell, i.e. — hsd + 
Hu, and by the so-called hybridization function A. 3 d{^) 
which describes the dynamic hybridization of the Co 3d- 
shell with the rest of the system (i.e. the bath). The 
hybridization function can be extracted from the non- 
interacting Green's function [G^] _1 = [Gc]^ 1 + ^3d(w) 
of the central region C as follows: 

A 3d M =LJ + fi-h 3d - [g^Hr 1 (4) 

gg d is the non-interacting Green's function projected onto 
the 3d-shell of the Co atom, i.e. gg d = PadGpPad. 
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Solving the generalized Anderson impurity problem is 
a difficult task, and at present there is no universal im- 
purity solver that works efficiently and accurately in all 
parameter regimes. Here we make use of impurity solvers 
based on an expansion in the hybridization strength given 
by A3d(aj) around the atomic limit. The starting point is 
an exact diagonalization of the (isolated) impurity sub- 
space i.e. the Co 3d-shell in our case given by the inter- 
acting Hamiltonian T-L^d- The hybridization of the impu- 
rity subspace with the rest of the system (given by the 
hybridization function A 3( ;(w)) is then treated perturba- 
tively. 

The so-called Non-Crossing Approximation (NCA) is a 
self-consistent perturbation expansion to lowest order in 
the hybridization strength. NCA only takes into account 
the most simple diagrams describing simple hopping pro- 
cesses where an electron or hole hops into the impurity at 
some time and then out at a later time (see Fig. [TJin App. 
|B|) . The so-called One-Crossing Approximation (OCA) 
improves upon the NCA by taking into account second 
order diagrams where two additional electrons (holes) are 
acommodated on the impurity at the same time in ad- 
dition to the NCA diagrams as shown in Fig. [7J OCA 
is the lowest order self-consistent approximation that is 
exact up to first order in the hybridization A$d oc V 2 
for all physical quantities while NCA is only exact to 
first order in A^d for the conduction electron self-energy. 
OCA improves considerably many of the shortcomings of 
NCA: It substantially improves the width of the Kondo 
peak and hence the Kondo temperature which now are 
only slightly underestimated. It also corrects the asym- 
metry of the Kondo peak. For very low temperatures 
(T <C Tjc), however, the height of the Kondo peak is still 
overestimated, and the Fermi liquid behaviour at zero 
temperature is not recovered. 

Hence, OCA is a reasonable approximation for solving 
the generalized impurity problem as long as the tempera- 
tures are not too low (i.e. more than one order of magni- 
tude below Tr-). In App. |B]we give a brief introduction 
to the NCA and OCA impurity solvers. A detailed de- 
scription of the NCA and OCA methods can be found 
e.g. in Refs. HSSM 



III. DFT ATOMIC AND ELECTRONIC 
STRUCTURE CALCULATIONS 

First, we have calculated the adsorption curves for a 
Co atom on a Graphene sheet with density-functional cal- 
culations within the generalized gradient approximation 
(GGA) with the CRYSTAL06H a b initio electronic struc- 
ture program for periodic systems employing an elabo- 
rate all-electron Gaussian basis set (6-31g). 

As possible adsorption sites for the Co impurity we 
consider the three high symmetry sites shown on the inset 
of the top panel of Fig. [TJ i.e. the hollow, top, and bridge 
site. For the calculations we have chosen a 5x5 unit cell in 
order to avoid interactions between impurities in different 
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FIG. 1: (color online) Adsorption of Co atoms on graphene. 

(a) Adsorption energy of Co adatom on graphene sheet in 
dependence of distance h for the three different adsorption 
sites as shown in the inset: Hollow (h), top (t), and bridge site 

(b) . (b) Magnetic moment of Co adatom on graphene sheet 
in dependence of distance h between Co atom and graphene 
sheet. 



unit cells. We have checked that the results essentially 
do not change when choosing a slightly smaller (4x4) or 
bigger (6x6) unit cell. Fig. []Ja) shows the adsorption 
energy for a Co atom in dependence of the distance h 
between the Co atom from the graphene sheet for the 
three different sites. One can see that the hollow site is 
the most favourable absorption site. We find the energy 
minimum for this site at a distance of h =I.5fA. As can 
be seen the energy minima for the other two adsorption 
sites are much more shallow. We find h — I.8lA for the 
bridge site, and h = I.96A for the top site. The top site 
is somewhat more favourable in energy than the bridge 
site. 

In Fig. [ljb) the magnetic moment of the Co atom 
in dependence of the height h of the Co atom over the 
graphene sheet is shown. Initially the magnetic moment 
is 3, i.e. that of a free Co atom. When the Co atom ap- 
proaches the equilibrium distance for a certain adsorp- 
tion site the magnetic moment starts to decrease. At 
the equilibrium distance, the magnetic moment of the 
Co atom has decreased to below 1.5 for the hollow site 
while for the bridge site it is more close to 1 already. 
For the top site , on the other hand, the magnetic mo- 
ment is still around 3 at the equilibrium distance as for 
the free atom. The magnetic moment for all three ad- 
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Ayco (e) 


fJ-Co 


A 1 




M3d 


h 


1.51 


-1.395 


0.536 


1.336 


1.127 


7.530 


1.197 


b 


1.81 


-0.660 


0.354 


1.175 


1.109 


7.533 


1.627 


t 


1.96 


-0.628 


0.248 


2.905 


2.913 


7.365 


2.150 


b* 


1.87 




0.311 


3.065 


2.890 


7.337 


2.150 



TABLE I: Summary of GGA results. Optimal height (h) 
of Co atom, adsorption energy (-E a d), charge (AQco) and 
magnetic moment (fJ.Co) of Co atom, total magnetic moment 
of unit cell (n), occupation (N^d) and magnetic moment fi^d of 
Co 3d-shell for the different adsorption sites considered here, 
i.e. h (hollow), b (bridge), t (top). In the last row (h*), we 
show the GGA results for the Co atom at the hollow site for 
the optimal distance h = 1.87A calculated with the B3LYP 
hybrid functional (see text). 



sorption sites approaches 1 when further decreasing the 
height of the Co atom. In the case of the bridge adsorp- 
tion site the transition from the free magnetic moment 
of 3 fiB to a magnetic moment of 1 /ib is remarkably 
abrupt. A summary of our GGA results is given in Tab. 
HI Our GGA results are in good agreement with previ- 
ous work by other groups using different method o 24 ' 39 . 
In particular, we find for the hollow adsorption site at 
the equilibrium distance that the magnetic moment of 
the 3<i-shell is indeed close enough to 1, implying a spin- 
1/2 for the Co 3c?-shell, as reported by Wehling et al^. 
However, the occupation of the 3c?-shell is around 7.5 and 
not around 9 as reported by Wehling et at As will be 
explained below the approximate spin- 1/2 of the Co 3d 
is essentially due to a hole in the Ei minority-spin levels. 

We have also performed hybrid functional 
calculations 41 which in general yield better geome- 
tries for molecules. Using the GAUSSIANO!) 4 ^ quantum 
chemistry code we have opimized the distance h of 
the Co atom on the hollow site of a finite graphene 
flake of hexagonal shape with 24 C atoms and with 
H atoms at the borders to saturate the bonds^ 2 -. Em- 
ploying the popular B3LYP hybrid functional together 
with the 6-31G basis set (also employed before in the 
CRYSTAL06 calculations) we now find a considerably 
higher value for the optimal height of the Co atom 
at the hollow site, namely h = 1.87A. For the GGA 
functionial we obtain very similar results (h — 1.47A) as 
before with CRYSTAL06 code. At the optimal height 
calculated with B3LYP (h = 1.87A) both, the GGA 
and the B3LYP results are remarkably similar with 
respect to the magnetic moments and occupations of the 
Co 3c?-shell: The occupation of the 3c?-shell is 7.4 and 
the corresponding magnetic moment is close to 2 now, 
implying a total spin of 1 for the 3<i-shell. In the last 
row of Tab. HI we have added the GGA results for the 
hollow site with the Co atom at the optimal distance 
h = 1.87A calculated with the B3LYP functional. 

In Fig. [2] we show the Co 3d-spectral functions calcu- 
lated with spin-polarized GGA for the Co atom at the 
hollow site for two different heights, i.e. the optimal 
height predicted by GGA (h — 1.51 A) and the optimal 
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(a) h=1.51A 




(b) h=1.87A 
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FIG. 2: (color online) Co 3d-spectral densities calculated with 
spin-polarized GGA for Co atom at hollow site for (a) the 
GGA equilibrium distance h = 1.51 A and (b) the B3LYP 
equlibrium distance h* = 1.87 A. Full (dashed) lines cor- 
respond to spin-up (down) states. The vertical dashed line 
indicates the Fermi level. 



height (h = 1.87 A) predicted by B3LYP. Due to the C6- 
symmetry of the system the Co 3d-shell splits into three 
symmetry groups: The non-degenerate Ai group consist- 
ing of the rotationally invariant 3d 3z 2_ r 2 -orbital only, the 
doubly degenerate Ei group consisting of the 3d xz - and 
3d yz -orbitals and the E2 group consisting of the 3d xy - and 
3d a .2_ 2/ 2-orbitals. The postitions of the peaks in the spec- 
tral functions show the effective energy levels for each of 
the symmetry groups. We can see that at the GGA equi- 
librium height (Fig. Ufa)) the spin- 1/2 is essentially due 
to a hole in the Ei -levels. The Ei minority-spin levels 
are exactly at the Fermi level and hence are half-filled, 
while the Ei majority-spin and the Ai and E2 minority- 
and majority-spin levels are well below the Fermi level. 
This level scheme seems to suggest a 3d 9 configuration 
as reported by Wehling et at 24 , although the true occu- 
pation of the 3d-shell is really only 7.5. The reason for 
this apparant discrepancy is the dynamic hybridization 
of the Co 3d-levels with the graphene substrate. 

The Co 3d-spectral density calculated at the B3LYP 
equilibrium distance is shown in Fig. [2f b) . The main 
difference with the spectral density at the GGA equi- 
librium (Fig. |2fa)) is that the peak corresponding to 
the minority-spin Ei -levels are now completely above the 
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FIG. 3: (color online) Imaginary part of the hybridization 
functions for the Co atom at the hollow site as shown in the 
inset for h — 1.87 A and for the three symmetry groups. The 
inset shows a graphene sheet with Co adatom at hollow site 
The red hexagon indicates the central region C containing 54 
carbon atoms. 



Fermi level, and hence the minority-spin Ei-levels are 
completely empty. Hence the 3<i-shell now yields a mag- 
netic moment of about 2, i.e. a spin-1 due to two holes 
in the Ei-levels. Note that now the total magnetic mo- 
ment of the Co atom is 3 due to a single electron in the 
Co 4s-orbital which is fully spin-polarized. Also note 
that in this geometry, both the Ai- and the 4s-level do 
not couple at all to the graphene substrate for symmetry 
reasons-^.. Our GGA results at the B3LYP equilibrium 
distance show some similarity to the GGA+U results of 
Wehling et alM 



IV. GGA+OCA CALCULATIONS FOR CO 
ATOM ON HOLLOW SITE 

The inset of Fig. [3] shows the central region C con- 
taining the Co adatom at the hollow site embeded into 
a perfect graphene sheet. The cell has been taken large 
enough in order to avoid interaction between Co atoms in 
the CRYSTAL periodic supercell calculation. The shape 
of the unit cell has been chosen to reflect the 6-fold sym- 
metry of the system. The C region can now be embeded 
directly into the graphene lattice set up by supercells of 
the same shape as C as described before in Sec. [TTJ 

Due to the C6-symmetry of the system the Co 3d- 
shell splits into three symmetry groups: The non- 
degenerate Ai group consisting of the rotationally in- 
variant 3d 3z 2_ r 2 -orbital only, the doubly degenerate Ei 
group consisting of the id xz - and 3(i az -orbitals and the 
E2 group consisting of the 3d xy - and 3d 2 ,2_ J( 2-orbitals. 

The right panel of Fig. [3] shows the hybridization func- 
tions for each symmetry group. As can be seen at low 
energies the 3G?3 Z 2_ r 2 -orbital (Ai) has an essentially flat 
and zero hybridization with the graphene sheet. This 





ef! (eV) 


Re A 3d (eV) 


Zf! (eV) 


Ax 


-3.650 


-0.144 


-3.794 


Ei 


-3.637 


+0.197 


-3.440 


E2 


-3.414 


-0.389 


-3.803 



TABLE II: Crystal field splittings for 3d-shell of Co atom 
on hollow site of graphene (h = 1.87 A). e 3( f are the on- 
site energies before double-counting correction, Re A^d is the 
real part of the hybridization function at the Fermi level, and 
= e^f + ReA3<j(0) are the effective on-site energies tak- 
ing into account the real part of the hybridization function. 
Note that the Ai and E2 effective energy levels are almost 
degenerate. 



can be understood by the fact that the coupling of the 
3d3 Z 2„ r 2 on the hollow site to the 2p z -orbitais of the 
graphene sheet is zero due to symmetry considerations^. 
The E2-group has the (overall) strongest coupling to the 
C 2p 2 -orbitals in that geometry as can be seen from the 
corresponding hybridization function. Note that the hy- 
bridization function for both the Ei and E2 symmetry 
group vanish linearly near the Fermi energy reflecting 
the linear DOS of the graphene sheet. Interestingly, the 
slopes for the Ei- and E2-hybridization functions are dif- 
ferent for energies below and above the Fermi level. Mo- 
rover, for negative energies the Ei-hybridization function 
has a steeper slope than the E2-hybridization, while for 
positive energies it is the other way round. This implies 
that Ei couples predominantly to the negative energy 
band of graphene while E 2 couples predominantly to the 
positive energy band. 

In Table. lU we show the energy levels, the real part 
of the hybridization function at the Fermi level, and the 
effective energy levels of the Co 3ci-shell. While the split- 
ting of the bare energy levels e%d is around 0.2 eV, the 
splitting of the effective energy levels e 3 d taking into 
account the (real part of the) hybridization with the 
graphene is somewhat bigger with around 0.4 eV. Taking 
into account the real part of the hybridization changes 
also the ordering of the energy levels. For the ordering 
of the effective energy levels we obtain: e(i?2) < e(^4i) < 
e(i?i). Note that this level ordering is similar to the one 
obtained by Wehling et a/.—. 



A. GGA+OCA results 

We now solve the generalized impurity problem of the 
Co 3<i-shell coupled to the graphene sheet by using the 
OCA impurity solver as described in Sec. Inland in App. 
|B1 We take the Coulomb repulsion U between <i-orbitals 
to be somewhat higher than in bulk due to the supposedly 
smaller screening: U — 5 eV. The Hund's rule coupling 
is taken to be the same as in bulk, i.e. J = 0.9 eV since 
it is less affected by screening effects. The DCC was 
calculated by the standard expression given by eq. ([3]) 
and taking the GGA value of about 7.5 for N34. With 
OCA, the total occupation of the impurity 3c?-shell is 
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FIG. 4: (color online) GGA+OCA results for U — 5 eV and J = 0.9 eV and total occupation N^d of the 3d-shell around 7.0. 
(a-c) Co 3d-spectral functions near the Fermi level for all three symmetry groups at different temperatures. The insets show 
the corresponding orbital occupations as a function of temperature, (d) Spectral function of p z -orbital of carbon atoms near 
the Co-atom for different temperatures. The black dashed line shows the p z -spectral function when the self-energy of the Co 
3d-electrons Y,$d is set to zero, (e) Comparison of GGA+OCA 3d-spectral density with spin-polarized GGA spectral density 
on a large energy scale. The inset shows a blow-up of the region near the Fermi level. 



now between 7.1 and 7.0 at low temperatures compared 
to the GGA value of 7.5. This result is very typical for 
correlated calculations: The correlations push the system 
towards integer occupation numbers to lower the energy. 

Fig. 2] shows the spectral functions of the Co 
3d-electrons for different temperatures calculated with 
OCA. At high temperatures we find peaks near the Fermi 
level for all three symmetry groups. However, the tem- 
perature dependencies of the individual peaks are quite 
different from each other: For the Ai-symmetry consist- 
ing of the d 3z 2_ r 2-orbital the peak becomes smaller at 
increasing temperatures and finally develops a step-like 
behaviour at low temperatures. The peaks in the Ei and 
E2-spectra on the other hand become more pronounced 
with decreasing temperature more or less as expected for 
Kondo effect. 

The form of the spectral functions and the temper- 
ature behaviour can be understood by looking at the 
occupation numbers for the different symmetry groups. 
The insets of Figs. HJa-c) show the occupation num- 
bers as a function of temperature for each of the three 
symmetry groups. For the Ai-symmetry the occupa- 
tion changes from sal. 6 at high temperatures to almost 
2 at low temperatures. Hence, at low temperatures the 
ti322_ r 2-orbital is practically completely filled, and exci- 
tations near the Fermi level are suppressed. For the Ei- 
symmetry group composed of the d xz - and c^-orbitals 
the occupation changes from 2.8 at high temperatures to 
2.1 at low temperatures. Hence at low temperatures this 
doubly-degenerate channel is essentially half-filled, cor- 
responding to a spin of almost 1 in this channel due to 
Hund's rule coupling. The peak near the Fermi level in 



the Ei-spectral function is hence due to a spin-1 Kondo 
effect in this channel. For the E2-symmetry composed 
of the d xy - and d a; 2_j / 2-orbitals we find an occupation 
of 2.8 for high temperatures, and of almost 3 for low 
temperatures. Hence the doubly-degenerate E2-channcl 
accomodates 3 electrons with a total spin- 1/2 so that 
the peak in the E2-spectral function is due to a spin- 
1/2 Kondo effect. Accordingly, the peaks near the Fermi 
level in the 3ci-spectral function at low temperature are 
due to a spin-3/2 Kondo effect in the 3c?-levels. This in- 
terpretation is further corroborated by the fact that the 
main contribution to the electronic configuration of the 
3rf-shell is an eightfold degenerate state with Nd = 7 and 
S = 3/2. The very fact that such a high spin state can be 
screened by the graphene substrate is an evidence of the 
multi-channel character of the graphene as a conduction 
electron bath. 

In Fig. 2Je) we compare the 3c?-spectral density calcu- 
lated with the GGA+OCA (at T = 120 K) method on 
the one hand and with spin-polarized GGA on the other 
hand. As can be seen the spectral densities calculated by 
the two methods differ significantly from each other. In 
the GGA spectrum, the peak just below the Fermi level 
(see also the inset of Fig. 0Je)) corresponds to a com- 
pletely filled i?i-shell with down-spin as can be seen by 
comparison with Fig. ^c). The peak above (and also 
a little farther from) the Fermi level corresponds to a 
completely empty spin-down i?2-shell. Both peaks carry 
essentially the full spectral weight of the correpsonding 
orbitals. In contrast, the two peaks near the Fermi level 
in the GGA+OCA spectra which are also of E\- and E2- 
symmetry, do not carry the full spectral weight of the 
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corresponding orbitals due to the renormalization of the 
quasi-particle by the electron-electron interactions which 
cannot be captured by a static mean-held calculation like 
GGA. 

Fig. HJd) shows that the correlations in the Co 3d- 
shell also affect the projected DOS of the p z -orbitals on 
the carbon atoms that are directly coupled to the Co 
atom (i.e. the 6 carbon atoms nearest to the Co atom). 
This can be understood as follows: Due to the coupling 
between the carbon p z - and the Co 3d-orbitals, the full 
p^-Green's function G v {ui) is related to the full Green's 
of the 3<i-orbitals Gd{oj) by: 

G p {u>) = G» + ]T Gl{uj)V pd G d {u:)V dp G p {uj) (5) 
d 

where G p (lu) is the Green's function of the (unperturbed) 
graphene p z -orbitals without the coupling to the Co 3d- 
orbitals and V p d = V^ p is the coupling between the p z - 
and individual 3<i-orbitals. The full 3<i-Green's func- 
tion Gd(uj) takes into account the correlations given by 
£3^(0;) and the hybridization with the p z -orbitals given 
by A3d(w), i.e. V p dGd(uj)Vd p represents the full T-matrix 
for the p z -electrons. 

One can see a temperature-dependent peak roughly 
at the position of the peaks in the Ei- and E2-spectral 
density. For comparison we also show the p z -DOS when 
the self-energy of the Co 3d-electrons £3,2 is set to zero. 
In this case the 3d-levels are far below the Fermi level 
and hence cannot affect the p z -DOS near the Fermi level. 
Hence the sharp resonance in the p z -DOS vanishes. This 
shows that the Kondo peaks in the Co 3<i-spectra can 
also be detected indirectly by spectroscopy of the nearby 
carbon atoms. 



B. Dependence on 3d-energy levels 

The double-counting correction term hd cc given by eq. 
([3]) which substracts out the Coulomb interaction within 
the correlated 3<i-subspace that has already been taken 
into account on a static mean-field level in the LDA or 
GGA calculation is not known exactly. Hence the ex- 
act position of the 3(i-energy levels is not known. Here 
we explore how shifting the 3<i-energy levels changes the 
results. 

Figs. OJa-c) show the Co 3<i-spectral densities for each 
of the different symmetries and for different shifts Ae^ of 
the 3<i-levels with respect to their original position, and 
Figs. [SJi,e the corresponding occupancies. Quite obvi- 
ously there is a drastic change in the correlations when 
the impurity-levels are lowered by more than 0.2 eV. At 
that point, a substantial peak appears in the Ai-spectra, 
while the peak in the Ei-spectra disappears. The peak 
in the E2-spectra on the other hand becomes qualita- 
tively different for energy shifts Ae^ < —0.2 eV, i.e. their 
spectral weight increases considerably. The changes in 
the spectra are accompanied by changes in the occupan- 
cies: Most importantly, the occupancy of the E2-levels 
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FIG. 5: (color online) 3d-energy level dependence, (a-c) Spec- 
tral functions of Co 3d-levels for different energy shifts Aa of 
the 3d-levels (as indicated in the legend of (c)) at T — 37 K 
and for the three different symmetries. (d,e) Orbital occu- 
pancies of the Co 3d-levels as a function of the energy shift 
for each of the three symmetries and for T — 37 K. 



increases abruptly by about 0.6 for Ae^ « —0.2 eV. 
The changes in the 3c?-levels with Ai- and Ei -symmetry 
are less pronounced but also occur at the same point. 
Also the overall occupation of the impurity levels changes 
quite abruptly by about 0.5 electrons from ~ 7 to 
Nd Ri 7.5. This indicates a strong change in the elec- 
tronic configuration of the Co 3d-shell. 

Indeed we find that for an energy shift around -0.2 eV, 
another 3c?-shell state with Nd — 8 and S = 1 becomes 
similar in energy to the state with Nd = 7 and S = 3/2. 
Hence, at this point the system enters the mixed-valence 
regime. This explains the strong increase of the spectral 
weight in the E2-peak for Ae<j < —0.2 eV since the mixed- 
valence regime is characterized by a strong peak near the 
Fermi level. On the other hand, the integer occupation 
of 2 in the Ei-levels, and the associated spin-1, essen- 
tially does not change upon entering the mixed-valence 
regime. In fact, the occupation is now even closer to 2 
than it was in the Kondo regime. Also the occupation 
of the Aj-level remains essentially constant upon enter- 
ing the mixed-valence regime. Hence the mixed-valence 
is actually only induced in the E2-levels since the occu- 
pations of the Ai-level and the Ei-levels in the Nd = 8 
state are the same as in the Nd = 7 state. The vanishing 
of the peak in the Ei-spectra upon entering the mixed- 
valence regime is due to the weaker hybridization with 
the graphene for negative energies (compare Fig. [3]): As 
the Kondo peak in the Ei-levels moves beyond the Dirac 
point as the impurity-levels are lowered in energy, the 
Kondo coupling in the Ei-levels becomes weaker and the 
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Kondo temperature much smaller so that the Kondo ef- 
fect is not observed anymore at the temperatures consid- 
ered here. It is not quite clear what triggers the peak in 
the Ai-spectra upon entering the mixed-valence regime. 
One possibility is that the increased charge fluctuations 
in the mixed valence regime also lead to charge fluctu- 
ations in this level via the Hund's rule coupling to the 
other 3<i-levels. The slight decrease in the occupation of 
the Ai-level in the mixed-valence regime supports this 
view. 

Note that the position of a peak (if present) changes ac- 
cording to the energy shift. This also explains the change 
in the width of the peaks for the Ei- and E2- symmetries: 
As a peak moves away from the Fermi level the broad- 
ening increases due to the increased hybridization with 
the graphene substrate. The Ai-level on the other hand 
does not couple at all to the graphene substrate (zero 
hybridization) so that the changes in the shape of the 
Ai-peak are purely correlation effects. 

Hence tuning the energies of the Co 3d-levels changes 
the correlations in these levels due to changes in the oc- 
cupancy. More precisely, the system is driven from the 
Kondo regime to the mixed- valence regime when the en- 
ergies of the 3<i-levels are lowered by more than 0.2 eV 
with respect to their original energies. The change in the 
correlations is reflected by corresponding changes in the 
3d-spectral functions. 



C. Gate voltage dependence 

In the previous subsection we have seen that chang- 
ing the exact position of the Co 3<i-levels with respect 
to the Fermi level has a considerable effect on the elec- 
tronic correlations and hence on the spectral density of 
the 3<i-levels. Although changing the impurity levels di- 
rectly is not feasible in an experiment, one can instead 
easily tune the chemical potential of graphene by appli- 
cation of a gate voltage. The chemical potential then is 
shifted with respect to both, the impurity levels and the 
graphene bands. Hence it should be possible to control 
the electronic correlations and in particular the Kondo 
effect by application of a gate voltage. The change in 
correlations can then be detected in the 3<i-spectra for 
example by STM spectroscopy. 

More precisely, a gate voltage v g shifts the chemical 
potential /i with respect to the graphene bands and the 
impurity levels according to (X = ej? — ev g where e denotes 
the elementary charge and cf is the Fermi level of neutral 
graphene, i.e. the energy of the Dirac points (see inset 
of Fig. [6^). Overall the results for shifting the chemical 
potential presented in Fig. [5] look very similar to the 
previous ones for shifting the impurity levels in Fig. [SJ 

At the hollow site only the Co 3<i-levels of Ei- and E2- 
symmetry couple to the graphene sheet for low energies 
while the Ai -level is completely decoupled. Hence only 
the spectra of Ei- and E2-levels can be measured by an 
STM. Therefore in Fig. H{a), we show the total spectral 



density of the Co Ei- and E2-levels for different values of 
the gate voltage. For gate voltages smaller than -0.1 eV 
there is only one peak in the spectral function which is 
due to the E2-levels. For gate voltages > —0.1 eV a sec- 
ond peak appears which originates from the Ei -levels. 
The E2-peak and also the Ei-peak (when present) move 
with respect to the chemical potential \i according to the 
change in gate voltage. Hence the peaks in the spectral 
function are not pinned to the chemical potential as ex- 
pected for normal Kondo effect, but are rather pinned to 
the Dirac points of the graphene sheet. For complete- 
ness we also show the spectral density of the Ai -level in 
Fig. [6p for different gate voltages. A peak in the spectral 
density at low temperature is only observed for negative 
bias voltages. The peak disappears for zero and positive 
gate voltage. Also the resonance in the Ai-level appears 
to be pinned to the Dirac point of graphene rather than 
to the chemical potential. 

As can be seen in Figs. [BJi,e also the individual occu- 
pations of the 3c?-levels as a function of the applied gate 
voltage show a similar behaviour as before as a function 
of the energy shift Ae^ of the 3d- levels. In fact, the expla- 
nation for the behaviour of the spectra and occupancies 
as a function of the energy shift of the 3d-levels also ap- 
plies to their behaviour as a function of the gate voltage. 
Here, when the gate voltage becomes smaller than -0.1 V 
the system leaves the Kondo regime with an Na = 7 
and S = 3/2 state as the principal contribution to the 
electronic configuration of the 3rf-shell, and enters the 
mixed- valence regime with an equal contribution coming 
from an Nd = 8 and S = 1 state. 

We can estimate a Kondo temperature from the width 
of the resonance in the Kondo regime. Fig. [5b shows the 
Kondo temperature of the E2-levels estimated from the 
resonance width T using the relation ksTn = n ■ w ■ T/8 
known from the Fermi liquid theory of the Anderson 
model where w is Wilson's number 0.4128 (see e.g. the 
book by Hewson^, chap. 5). In Fig. [U we have also in- 
cluded "Kondo temperatures" estimated from the width 
of the peaks for negative gate voltages where the system 
is not in the Kondo regime anymore. Hence outside the 
Kondo regime it is only a measure for the width of the 
resonance near the chemical potential. 

Note that in the Kondo regime (i.e. for v g > 0) the 
width of the resonance and hence the Kondo temperature 
are essentially constant. This behaviour is different from 
the one before where we only had shifted the impurity 
levels. There the width of the Ei- and E 2 -peaks increased 
as they moved away from the Fermi level due to the linear 
increase of the hybridization with the graphene substrate. 
Here, however, the peak is pinned to the Dirac point of 
graphene, i.e. it moves together with the graphene bands 
with respect to the chemical potential. Hence the hy- 
bridization with the graphene conduction electrons does 
not change here. The strong dependence of the width 
of the resonance on the gate voltage for negative gate 
voltages on the other hand is due to the changes in the 
electronic correlations as the system is driven out of the 
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FIG. 6: (color online) Gate voltage dependence, (a) 3d-spectral function for Ei and E2 symmetry for different gate voltages 
v g (T = 37 K). Different offsets have been added to the spectral functions of different v g in order to better distinguish them 
from each other. The inset shows a Dirac cone of the graphene sheet with the position of the chemical potential jj, for a certain 
gate voltage v g . (b) Spectral function for Co 3d-level of Ai-symmetry for different gate voltages v g (T — 37 K). The chemical 
potential in both (a) and (b) is always at zero frequency, (c) Kondo temperature estimated from width of E2 resonance as 
a function of gate voltage v g at T = 37 K. (d-e) Occupancies of Co 3d-levels resolved for the three different symmetries as a 
function of gate voltage v g at T = 37 K. 



Kondo regime. 



V. SUMMARY AND CONCLUSIONS 

Using DFT electronic structure calculations we have 
found that a single Co atom most likely adsorbs at the 
hollow site of an otherwise clean graphene sheet. Fur- 
thermore the DFT calculations show, that the magnetic 
moment and the electronic configuration of the 3d-shell 
of the Co atom depends quite strongly on the adsorption 
height over the graphene sheet. At large distantes the 
formation of a local moment is favoured while at short 
distances the magnetic moment is quenched due to strong 
hybridization with the substrate. Since the adsorption 
energy curves are quite shallow, it might be possible in 
an experiment to push the Co atom from the local mo- 
ment to the weakly correlated regime. 

Using our GGA+OCA correlated electronic structure 
method for nanoscopic systems, we have then studied 
the impact of dynamical correlations due to the strongly 
interacting Co 3d-clcctrons on the electronic structure. 
Because of the strong crystal field splitting of the 3d- 
levels by graphene substrate, the behaviour of each of the 
Co 3d- levels strongly depends on its orbital symmetry. 

We find that the dynamic correlations can give rise 
to a spin-3/2 Kondo effect. The spin-3/2 is composed 
of a spin-1 due to two electrons within the doubly- 
degenerate orbitals of Ei-symmetry and a spin- 1/2 due 
to a hole within the doubly-degenerate orbitals of the E2- 



symmetry. The Kondo effect gives rise to corresponding 
peaks near the Fermi level in the Ei- and E2-spectral 
function. The Ai-channel on the other hand is nearly 
completely filled and hence in the Kondo regime, it does 
not give rise to a peak near the Fermi level. The real- 
ization of a spin-3/2 Kondo effect points to the multi- 
channel character of graphene. Hence a spin- 1/2 impu- 
rity could in principle give rise to multi-channel Kondo 
physics and Non-Fermi liquid behaviour. 

The Kondo effect can be controlled by a gate voltage 
which changes the chemical potential and thereby the 
filling of the 3d-levels. The change in the correlations 
can be observed by a stark change in the spectra as the 
system is driven out of the Kondo regime and into the 
mixed valence regime. Another remarkable result is that 
the Kondo peaks appear to be pinned to the Dirac point 
of the graphene substrate rather than to the chemical po- 
tential when a gate voltage is applied. This last finding is 
actually in agreement with very recent STM spectroscopy 
measurements of Co atoms on graphene which find peaks 
in the dl/dV curves that move according to the applied 
gate voltage^. But it has not yet been pointed out by 
other theoretical work considering the dependence of the 
Kondo effect in graphene on a gate-voltag o 15 i 26 . 
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Appendix A: Details of the embedding method 

Here we show how to compute the embedding self- 
energy Eh given in eq. ([lj that describes the dynamic 
hybridization of the central region C with the rest of the 
system (i.e. the graphene sheet). First, we obtain an ef- 
fective mean-field description of the host material from a 
GGA calculation for a supercell of clean graphene corre- 
sponding to the C region. From the Kohn-Sham Hamil- 
tonian Hh (k) of the host material in reciprocal space we 
calculate the fc-dependent Green's function: 

G B (k,u) = (u + f i-H B (k) + i v )- 1 (Al) 

The local KS Hamiltonian of the supercell is obtained by 
summing up Hn(fc) over all k-points within the Brillouin 
zone: 

H h ,o = ]THh(£) (A2) 

k 

Similarly, the local Green's function of the supercell is 
given by summing up Gn(fc,w) over all k-points within 
the Brillouin zone: 

Gh,o(w) = ^G H (fc,w) (A3) 
k 

On the other hand, we can write the local Green's func- 
tion also in terms of the local Hamiltonian Hh.o of the 
supercell and the embedding self-energy (which describes 
the dynamic hyrbidization with the rest of the graphene) 
as: Gh,o(w) = + ^ — Hh,o — Sh(w)] -1 . Hence the 
embedding self-energy is given by: 

EhM =w + /i-H H ,o - [Gh^MP 1 (A4) 

Appendix B: The NCA and OCA impurity solvers 

The general multi-orbital Anderson impurity model 
can be written in the following form: 

n = ]T e d a dld a + I U aMs dl$d n d 5 (Bl) 

a/3 afl-yS 

^kvC^u^ku 

kvota kva 

where in order to keep the notation simple we have com- 
bined the spin- and orbital degrees of freedom into one 
index for each impurity level a and each band v. 



The Non-Crossing Approximation (NCA) and the 
One-Crossing Approximation (OCA) both solve the An- 
derson impurity model by expansion in the hybridization 
strength around the atomic limit. The starting point is 
an exact diagonalization of the impurity subspace, i.e. of 
the Co 3d-shell in our case, including the Hubbard-like 
interaction term: 

^|m)£ m <m| (B2) 

m 

where |m) are the many-body eigenstates oiT-L^d and E m 
the respective eigen energies. 

One now introduces auxiliary fields a m , a' m (called 
pseudo-particles) such that each impurity state is rep- 
resented by a corresponding pseudo-particle: 

4JPPV) = |m) (B3) 

where |PPV) is the pseudo-particle vacuum. The com- 
pleteness of the impurity eigenstates imposes the follow- 
ing constraint: 

Q = Y,^ m = l (B4) 

ni 

The physical electron operators ct^ can now be expressed 
by the PP operators: 

dl = ]T(F Qt ) Ilm 4a m (B5) 

where (F a ^) nm = (n\dJ a \m) are the matrix elements of 
the impurity-electron creation operator. For later con- 
venience we also define the corresponding matrix ele- 
ments of the impurity-electron destruction operator as: 
(F a ) nm = (n\d a \m). The anti-commutation rules for the 
physical electron operators then require that the PP a m 
is a boson (fermion) if the corresponding state |m) con- 
tains an even (odd) number of electrons. 

In the PP representation we can now rewrite the 
Hamiltonian of the generalized Anderson impurity model 
as follows: 

ft = ^ E m al n a rn + e kv c\ v c kv + A(Q - 1) 

m kv 
+ J2{ V ^A V & l( Fa )nm&n + H.C.) (B6) 

kua 

where we have included the constraint Q = 1 into the 
Hamiltonian. The corresponding Lagrange multiplier A 
can be interpreted as a (negative) chemical potential for 
the PPs. 

In the PP picture, the hybridization with the bath elec- 
trons given by the last term in eq. (|B6j) becomes now the 
interaction for the PPs. Because of the fermionic and 



11 




FIG. 7: Diagrams for pseudo-particle self-energies in NCA 
and OCA approximations. 

bosonic commutation rules for the PPs, one can now de- 
velop a diagrammatic perturbation expansion in the PP 
interaction. The PP propagators can then be written as 

G m {u) = (u-X-E m -H m (u))- 1 (B7) 

where S m is the PP self-energy describing the dynamic 
interaction with the other PPs. 

The Non- Crossing Approximation consists in taken 
into account the diagrams shown in the first row of Fig. 
[7] for a certain pseudo-particle. The NCA diagrams de- 
scribe processes where a single electron (hole) jumps 
from the bath to the impurity and back thereby tem- 



porarily creating a PP with N+l (N-l) electrons. The 
NCA equations correspond to a self-consistent perturba- 
tion expansion to lowest order in the hybridization func- 
tion A a (w) = J2k v ^kv a^kv.a- Since the fermionic self- 
energies depend on the dressed bosonic propagators, and 
vice versa, the NCA equations have to be solved self- 
consistently. Once the NCA equations are solved the 
physical quantities can be calculated from the PP self- 
energies. 

OCA takes into account second order diagrams where 
two bath electron lines cross as shown in the right col- 
umn of Fig. [7J The self-energies for the fermions 
(bosons) again depend on the full propagators for bosons 
(fermions), and hence the OCA equations have to be 
solved self-consistently. The explicit expressions for the 
OCA self-energies are second in the bath hybridization 
function A a . The OCA is the lowest order self-consistent 
approximation that is exact up to first order in the hy- 
bridization A^d oc V 2 for all physical quantities while 
NCA is only exact to first order in A 3( ; in the PP and 
conduction electron self-energies. Further details of the 
NCA and OCA impurity solver can be found e.g. in Refs. 
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